COMBINATORICS OF LABELLED PARALLELOGRAM 

POLYOMINOES 

J.-C. AVAL, F. BERGERON, AND A. GARSIA 



Abstract. We obtain explicit formulas for the enumeration of labelled parallelogram 
polyominoes. These are the polyominoes that are bounded, above and below, by north- 
east lattice paths going from the origin to a point (k,n). The numbers from 1 and n (the 
labels) are bijectively attached to the n north steps of the above-bounding path, with 
the condition that they appear in increasing values along consecutive north steps. We 
calculate the Frobenius characteristic of the action of the symmetric group S„ on these 
labels. All these enumeration results are refined to take into account the area of these 
polyominoes. We make a connection between our enumeration results and the theory of 
operators for which the intergral Macdonald polynomials are joint eigenfunctions. We 
also explain how these same polyominoes can be used to explicitly construct a linear 
basis of a ring of £> .^-invariants. 
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1. Introduction 

Parallelogram polyominoes have been studied by many authors (see [5, 10, 8] for a nice 
survey and enumeration results). They correspond to pairs tc = (a,j3) of north-east 
paths going from the origin to a point (k, n) in the combinatorial plane NxN, with the 
path a staying "above" the path (3. Our aim here is to study properties, and related nice 
formulas, of "labelled parallelogram polyominoes". These are obtained by bijectively 
labelling each of the n north steps of the path a with the numbers between 1 and 
n. Our motivation stems from a similarity between this new notion and recent work 
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on labelled intervals in the Tamari lattice, in connection with the study of trivariate 
diagonal harmonic polynomials for the symmetric group (see [4]). 

We calculate explicitly the Frobenius characteristic of the natural action of the symmetric 
group S n on these labelled polyominoes; and study aspects of a weighted version of this 
Frobenius characteristic with respect to the area of the polyominoes. This connects our 
study to interesting operators for which adequately normalized Macdonald polynomials 
are joint eigenfunctions. This is the same theory that appears in the study of the S n - 
module of bi-variate diagonal harmonics (see [12, 13, 15]). 

We also extend some of our considerations to "doubly" labelled parallelogram polyomi- 
noes, with added labels on east steps of the below-bounding path; with a corresponding 
action of the group §i; x S„. Several components of these spaces are naturally related to 
parking function modules. 

2. Paths and labelled paths 

A k x n north-east (lattice) path in N x N is a sequence a = (p , . . . ,Pi, . . .p^) points 
Pi = (x{, Hi) in N x N, with p = (0, 0), such that 

{(xi, Hi) + (1.0) an east step, or 
\ x iiUi) + (0, 1) a north step. 

We denote by Pk,n the set of north-east paths from (0, 0) to (k, n) (i.e. pjq = (k, n), and 
hence N = k + n). As is well known, these paths number ( n ^ fe ). They are often encoded 
as word u = W\w2 ■ ■ -Wn, on the alphabet {x, y}, with x standing for an east step, and 
y for a north step; and k is the number of x's, while n is that of y's. Thus, the path of 
Figure 1 is encoded as yyxxxyyxxyxxxyxx . 

Another description of such a path a may be given in terms of the sequence of "heights" 
of its k horizontal steps (written from left to right), so that we may write 

a = a\a 2 ■ ■ ■ a^. 

More specifically, aj is equal to the unique yj such that Xj = % and Xj-i — % — 1. We say 
that this the height sequence description of a. Thus, the top path in Figure 1 may be 
encoded as the height sequence 2224455566. Perforce, the height sequence is increasing, 
i.e. dj < a i+ i, and any increasing sequence a\a 2 ■ ■ ■ a k , with < a, < n, corresponds to 
a unique path. 

Exchanging the role of the axes, we may also describe a in terms of its indentation 
sequence. This is simply the bottom to top sequence of distance between vertical steps 
in a, and the y-axis. For example, the indentation sequence of our running example is 
003358. Once again, this establishes a bijection between k x n north-east lattice paths, 
and length n increasing sequence of values between and k. 

There is a classical bijection between north-east paths in Pk,n an d monomials of degree 
k in the variables x = xi, . . . ,x n+ i. One simply associates to a height sequence a = 
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(di, a,2, ■ ■ ■ , dk), the monomial 

Thus, we get the well known formula for the Hilbert series, denoted R(x), of the 
polynomial ring R = C[x]: 



R(x) 



k=0 



1 — X 



n+1 



(1) 



(2) 



Recall that the coefficients of this series correspond to the dimension of the homogeneous 
components of the space considered. 



The q- weighted area enumeration of Pfc jn , is defined to be 

P*,n(9) := <? arCa(a) ' 



(3) 



with area(a) equal to the number of cells inNxN lying below the path a. In Figure 1, 
the area of the path is the shaded region. It is also well known that P& jn (?) is equal to 
the classical q- analog of the binomial coefficient: 

n + k 
k 



Pfc,n(g) 



J q 



h k (l,q,...,q n+1 ) 



(4) 
(5) 



n < 



"V - 

k 



Figure 1. A labelled path, and its area. 



Labelled paths. We now add labels to north steps of paths in Pfc )Tl . We require these 
labels, going from 1 to n, to be increasing whenever they lie on north steps having same 
horizontal coordinate. The underlying path, of a labelled path £, is said to be its shape. 
It is easy to see that the number of labelled paths is (k+l) n , with g-weighted enumeration 
given by the straightforward g-analog 

([*+l]g)"= E ( 6 ) 
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We denote by hk, n , the set of labelled paths having shapes lying in Pk, n - Observe 1 that 
labelled path i may be bijectively turned into a function from [n] := {1,2, ... ,n} to 
[k + 1] := {1, 2, . . . , k + 1}, as follows. We simply set E(i) := j, whenever i lies along a 
north step having first coordinate equal to j — 1. Thus, the fibers of the resulting function 
correspond to the labels of consecutive north steps. We say that the set partition of I, 
is the partition of [n] into these fibers. 



n < 



k 



Figure 2. A labelled path with below area equal to 41. 



The symmetric group S n acts by permutation of labels, up to reordering labels that 
lie one above the other. In fact, it may be best to consider this in the context of the 
indentation sequence encoding of paths. In this point of view, a labelled path is simply 
a permutation of this sequence. Thus, k x n labelled paths are just another name for 
length n sequences of numbers between and k. The shape a (underlying path), of a 
labelled path £, is simply the increasing ordering of this sequence. Still we will keep using 
the path terminology, which is better adapted to our upcoming study of parallelogram 
polyominoes. 

It is classical that the corresponding Frobenius characteristic is given by the formula 2 

Frob(L fc ,„)(z) = h n [(k + l)z] (7) 

= ^k + iy^EM ( 8 ) 

fi\-n ^ 

= ^T^(M,.^)s M (z). (9) 



fc+1 



Recall that this last equality says that 

s M (* + l):=s M (l, !,...,!) 



fc+i 



Our reason for considering functions from [n] to [k + 1] as labelled paths will become apparent in 
the sequel. 

2 We follow Macdonald's notation (see [16]) for symmetric functions. We also use plethystic notations. 
See [3] for more on this. 
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is the multiplicity, in Lk >n , of the irreducible representation indexed by the partition /x. 
Taking into account the area, we have a direct g-analog of the previous formula 



Frob(L fej „)(z; q) = h r 



Q 



fe+i 



1 



q-1 



E l n M+l)k^ (z) 

^ 2 M A c [l\q 

^8^(1, q,..., q k ) s M (z). 

fi\-n 



(10) 

(11) 

(12) 



3. Parallelogram polyominoes 



Beside setting up notations, the aim of this section is to recall basic (well-known) facts 
about parallelogram polyominoes. For two paths a and /3 in Pfc jTl , we get a polyomino 
7i = (a, j3) if a stays "above" /?. One usually thinks of a parallelogram polyomino 
tt = (a, (3) as the region of the plane bounded above by the path a, and below by the 
path (3. Except for endpoints, all points of the path f3 = (g , . . . , g i; . . . q N ) are required 
to be strictly below those of a = (j>o, . . . ,pi, . . .pn)- Hence, writing = (x^, y^), this is 
to say that y[ < yj whenever x\ = Xj, for 1 < i, j < N. We say that n is the height of 
tt, and that k is its width; and we denote by the set of parallelogram polyominoes 
of height n and width k. The number of 1 x 1 boxes lying between the two paths is 
always larger or equal to k + n — 1, and equality holds for ribbon shapes. Thus, it is 
natural to subtract this value in the definition of the area of the polyomino tt, denoted 
by area(7r). Hence ribbon shapes have zero area. The area-enumerating polynomial 
is then defined as 

P*,„(<?) := Yl <7 ar ° aW - ( 13 ) 

Clearly the reflection in the diagonal x = y maps bijectively the set Pfc jTl on the set P n> fe, 
we have the symmetry P^. n (g) = P n fc(g). This is reflected in the generating function 

F(x,y;q) := £ ¥ k>n (q) x k y n , 

(fc,n)GNxN 

= xy + xy + x y + xy + (2 + q)x y + x y 

+xy 4 + (3 + 2 q + q 2 )x 2 y 3 + (3 + 2 q + q 2 )x 3 y 2 + x A y 
+xy 5 + (4 + 3g + 2g 2 + q 3 )x 2 y 4 

+ (6 + 6 q + 5 q 2 + 2 g 3 + q 4 )x 3 y 3 

+ (4 + 3 q + 2 q 2 + q 3 )x 4 y 2 + x 5 y + . . . 
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which can be expressed in terms of Bessel functions (see [5]). When q — 1, we have 



det 



IPfc+l,n(l) — |Pfc+l,re| — s k,k( l] — 

n+1 

h k [n + l] h k+1 [n + l] 
hk-i[n + l] h k [n + l] 

n + k\ 2 fn + k + l\fn + k -1 

k ) V k + 1 )\ k-1 
1 /n + k\/n + k-l 
k+l{ k ){ k 



(14) 

(15) 

(16) 
(17) 



Indeed this may be readily deduced using Lindstrom-Gessel-Viennot [11] approach to 
the enumeration of non-intersecting path configurations in terms of determinants. 
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Figure 3. A Parallelogram polyomino and the corresponding Motzkin path. 



Alternate descriptions. We may shed new light, on our study of parallelogram poly- 
ominoes, by exploiting alternate encodings for them. For paths a, (3 G L fc+l n , respec- 
tively 3 given as increasing sequences a = a^ai • • ■ a k and = &162 • • ■ the pair 
7r = (a, (3) is a parallelogram polyomino if and only if ao = 0, b k +i — n+1, and a« > bi, 
for 1 < i < k. Hence, we may identify parallelogram polyominoes with semi-standard 
tableaux (see Figure 4) of shape k 2 , with values in{l,2,...,n + l}, thus explaining 
formula (14). Naturally, we may also describe parallelogram polyominoes in terms of 
indentation sequences, thus getting semi-standard tableaux of shape n 2 , with values in 
{1,2, ... ,n}. Going from the first encoding of a given polyomino to the second one, es- 
tablishes a classical bijection between the two equally numerous families of semi-standard 
tableaux. 



<2i 


a-2 


a 3 




ak 


h 


b 2 


h 




bk 



Figure 4. Semi-standard tableau of shape k 2 ; ai < aj + i and > bi. 



'The shift in indices is intentional here. 
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There is also a well-known bijection (see [8]) between parallelogram polyominoes and 
Primitive Motzkin Paths. Recall that these are the paths, in N x N, going from (0, 0) 
to (N, 0), with steps either north-east, east (either red or blue), or south-east. Thus, the 
consecutive points, along the path, are linked as follows 



(x i+ i,y i+ i) 



(xi, i/i) + (1, 1) a north-east step, or 

(xi, Ui) + (1, 0) a red-east step, or 

(xi, yi) + (1, 0) a bleu-east step, or 

{xi,yi) + (1, —1) a south-east step. 



Primitive paths are those that never return to the horizontal, except at endpoints. 
Motzkin paths are often simply presented in terms of Motzkin words. Recall that 
these are the words u = W1W2 ■ ■ ■ wn, on the alphabet {d, d,r, b}, such that 



(1) 
(21 



\ u <i\d - k<*ld' for all 1 < i < iV - 1, 



M 



Wi 



where \u\d (resp. |w<j|^) stands for the number of d's occurring in uj (resp. d), and u<i 
denotes the prefix of length i of u. The word corresponds to a primitive path, precisely 
when all inequalities are strict in the first condition. Each letter encodes one of the 
possible steps: d for north-east, r for red-east, b for blue-east, and d for south-east. For 
a polyomino 7r = (a, 0), with a = UiU 2 ■ ■ ■ and (3 = a = V\v 2 ■ ■ -v^ given as words in 
{x,y}, Viennot's encoding [8] consists in setting 

d if (ui,Vi) = (x,y), 
r if (ui,Vi) = (y,y), 
b if (ui,Vi) = {x,x), 
d if (ui,Vi) = (y,x). 
See Figure 3 for an example. 

Lastly, a parallelogram polyomino may be encoded as a word in the ordered alphabet 

A = {0,1,T,2,2,3,3,4,---}. 

To do this end, we start with its primitive Motzkin path encoding ir, and successively 
replace each step in tt by a m-sequence of letters in the A, < m < 2, according to the 
following rules. 

• A down step is replaced by an empty sequence; 

• an up step, at height j, is replaced by the 2-sequence j (j + 1); 

• a red horizontal step, at height j, is replaced by j; 

• a blue horizontal step at height j is replaced by j. 



Figure 5 illustrates this process. 
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FIGURE 5. From Motzkin path to word in A. 



Remark 1. The correspondence between parallelogram polyominoes and primitive Motz- 
kin paths, suggests that we should study "sequences" of parallelogram polyominoes which 
correspond to general Motzkin paths. In particular, this will make apparent more ties 
with the study of parking functions that correspond to Dyck path touching the diagonal 
at specific points. 



Figure 6. A parallelogram polyomino sequence. 



4. Parallelogram polyominoes, as indexing set of SL 2 -invariants 



One of our motivation, for our current study of the combinatorics of parallelogram poly- 
omino, is to tie this study to that of parking functions and labelled intervals in the 
Tamari poset. Recall that more and more evidence shows that this combinatorics is 
intimately related to the study of the § n -modules of bivariate and trivariate diagonal 
harmonics. Along these lines, it is interesting to observe that the number of polyomi- 
noes corresponds to the dimension of the space of SL 2 -invariants of weight k + 1 of the 
Grassmanian G(2,n + 1), see [18, page 238]. We make this more explicit as follows. Let 
X be the matrix of 2 x n variables 

\Vl, V2, ■■■ ,Vn) ' 

The special group acts on polynomials f{X), in C[X], by left multiplication of X by 
matrices in SL 2 . Invariants under this action are precisely the polynomial expressions 
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in the 2x2 minors 4 of X. To construct an explicit linear basis of the resulting ring, 
TZ n := C[X] SL2 , one needs only take into account the Pliicker relations on these 2x2 
minors. More precisely, as described in [17], we have the exact sequence 

— ► V n — ► C\Yi ,j\l<i<j<n — ► K n — ► 0, (18) 

where the Yij are variables, the middle arrow sends to JTjj (the 2x2 minor corre- 
sponding to the choice of columns i and j in X), and V n stands for the ideal generated 
by the Pliicker relations 

Y u Y jk -Y ik Y jt + Y ij Y ke = 0, for 1 < % < j < k < £ < n. (19) 

Moreover, using a somewhat reformulated result also discussed in [17, Thm 14.6], it fol- 
lows that a linear basis of lZ n = C[X] SL2 may be indexed by parallelogram polyominoes. 
In fact, we do this in a graded fashion 5 , constructing for each d > the ci-homogeneous 
part of the basis as a collection of product of d minors. The characterizing property 
of these products is that they must not be divisible by a binomial X^Xj^, for {i,j} 
incomparable to {k,£} relative to the following order on pairs of integers. We set 

{i, j} -< {k,£}, iff i < k and j < £, 

assuming, without loss of generality, that i < j and k < £. 

Now, given a polyomino n = (a,f3) 6 Pd+i,n-i, for which the height sequences of a and 
(3 are respectively 6 a^ai • • - aa and 6162 ■ • • bd+i, we consider the minor monomial 

d 

:= JJX( 6i+ i i(li+ i). (20) 

i=l 

The fact that the path /3 remains below the path a is equivalent to the fact that the 
minor monomial X n satisfies the charactering property of the previous paragraph. For 
example, the minor monomial associated to the polyomino ir = (2224455566, 0011112223) 
of Figure 7 is 

7v — A 13 A 23 A 25 A 36 A 47- 

Proposition 1. The family {X^^, with n varying in the set of parallelogram polyomi- 
noes Pd+i,n-i; constitute a linear basis for the degree d homogeneous component of the 
ring TZ n . 

It follows that, for a given n, the Hilbert series of TZ n is given by the formula 

oc 

TZ n (x) = J2\ F ^\x k (21) 



d=0 

n-2 



I _ x 2n-l k+ I I k )( k ) 

k=0 V 7 V 7 



x\ (22) 



4 These are clearly SL2-invariants. 
5 Considering 2x2 minors to be of degree 1. 
6 The shift in indices is intentional here. 
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This approach may readily be expanded to cover families of r non-intersecting paths from 
(0, 0) to (k, n), which also parametrize SL r -invariants. The corresponding enumeration is 
easily obtained in terms of Schur functions (using Lindstrom-Gessel-Viennot [11]), giving 



»(*■) 

fc + l,7ll 



(23) 

n+1 

where Pj^ denotes the set of such families, and k r stands for the rectangular partition 




Action of §„ on lZ n . We may refine the graded enumeration of the ring lZ n , by consid- 
ering its irreducible decomposition with respect to the action of S n , which corresponds 
to permutations of the n columns of X. Thus we consider the graded Frobenius charac- 
teristic 

oo 

Frob(ft rt )(z; x) = Frob (^i d) )( z ) x d , ( 24 ) 

d=0 

where Tlffl stands for the degree d homogeneous component of lZ n . 

To calculate this, we follow an approach due to Littlewood for the calculation of the 
character of the restriction to § n of an action of GL n . It follows from general principles 
that the character of the action of GL n on IZn , the degree d + 1 component of TZ n , 
is equal to the schur function s^(gi, . . . , g ra ). The corresponding character value, for 
the restriction of this action to S n , is obtained by evaluating this Schur function at the 
eigenvalues of permutations. This may be easily calculated, without explicit knowledge 
of these eigenvalues, in the following manner. 

We start by expanding Sdd(qi, . . .,q n ) in terms of the power-sum symmetric function, 
thus turning the required evaluation into the calculation of the sums of fc th -powers of the 
eigenvalues of any given permutation matrix a. But this turns out to be equal to the 
number of fixed points of the fc th -power of a. If \i is the partition giving cycle sizes of a, 
this number of fixed points may simply be expressed as 

fix(a fc ) = ^/i 4 . (25) 

fii\k 

Hence, the Frobenius characteristic that we are looking for is given by the formula 

n^ +1 \z) = y Sdd \ (26) 

where we set fk(^) '■= J2^\k^iy J '- e - the sum of the parts of \i that divide k. 

5. Labelled parallelogram polyominoes 



A labelled parallelogram polyomino is obtained by labeling (as in Section 2) the 
top path of a parallelogram polyomino. Equivalently, this may be presented in terms of 
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labelled Motzkin paths, with labels associated with either north-east steps, or red- 
east steps. In this case, the relevant increasing-labelling condition simply states that 
labels should be increasing whenever we they occur on consecutive north-east or red-east 
steps. Figure 3 gives an example of such a labelling, in both encodings. 
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Figure 7. A labelled parallelogram polyomino. 



We say that the underlying parallelogram polyomino is the shape of the labelled paral- 
lelogram polyomino, and denote by the set of labelled parallelogram polyominoes of 
shape 7r. We also denote by h k>n the set of labelled parallelogram polyominoes of height 
n and width k, so that 

A maximal length sequence of north steps, in a path a, is said to be a big-step of a. 
Clearly, for a height n path a, the sequence of big-step lengths in a forms a composition 
of n. We denote it by 7(a), also setting 7(7?) := 7(a) in the case of a parallelogram poly- 
omino 71 = (a, (3). Recall that it is customary to write u \= n when v is a composition 
of n. 

Remark 2. It is surely worth mentioning that, when k = n, "classical" Parking Func- 
tions may be obtained as special instances of parallelogram polyomino (sequences'), in 
fact in at least three different ways. The first is obtained by the simple device of con- 
straining the lower path to be equal to the zigzag path, the one who has height sequence 
equal to 012 • • • (n — 1). In a similar manner, we may fix the lower path to be the one that 
follows the lower boundary of the surrounding rectangle. Still, an even more natural way 
is to consider labelled parallelogram polyominoes sequences whose shape is symmetric 
with respect to the diagonal x = y. All of these correspondences are compatible with 
the action on labels, and area is easy to track. Hence, most of the following formulas 
may be considered as extensions of analogous formulas for parking functions, potentially 
in three different ways. 

For a height n parallelogram polyomino tt = (a, (3), the symmetric group S n acts by 
permutation on the set L^, up to a reordering of labels along big-steps. This implies 



This is needed to take into account the possibility that the two paths may touch one another. 
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that 



where we use the classical multinomial coefficient notation 
/ n\ n 



U\=[ ^^), (27) 



when v = (u t , u 2 , . . . , u k ). 

vj vi\h> 2 \ ■ ■ ■ v k \ 
It follows (see Proposition 3) that 

Proposition 2. The number of labelled parallelogram polyominoes is given by the formula 

l L *.»l=*' l " 1 ( n t-l 2 )' (28) 

and we have the identity 

It may be worth mentioning that we have the nice generating function 

n _x (n + k - 2\ ( 1 



n=0 



The following result provides a common refinement of formulas (17) and (28), stated in 
terms of the Frobenius characteristic of the S n -modules and "Lk,n (here considered as 
free vector spaces over the specified sets). 

Writing v = (u±, z/ 2 , . . . , Uj) for 7(71), it is classical that the Frobenius characteristic of 
is simply given by the symmetric function h u (z) (in the variables z = zi, z%, z%, . . .). 
Indeed, the Young subgroup x • • • x E> u . is precisely the stabilizer of the labelled 
parallelograms polyominoes of shape n. We thus have that 

Proposition 3. The Frobenius characteristic ofLk, n is given by the formula 

Frob(L fci „) := ^ Kw( z )' (31) 



and 



Frob(L M ) = ^h n [kz] y k k _ l J (32) 

j2 k m-i( n + k _- 2 )EM (33 ) 



Ahn 




(34) 



' V., 1 ! ^1 • • • 1 

k 

where l[y) (resp. £(fj)) stands for the length of a composition v {resp. partition fi), and 
k is considered as a constant for the plethystic evaluation in the formulas. 
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Observe that Formulas (28) and (17) follow easily from the above expressions, and 
that (32)-(34) are polynomial in k. Special values are as follows: 

Frob(L M ) = ^(z) 
'k\ 

)h n (z) + kh 2 (z) 



Frob(L fc , 2 ) 
Frob(L fci3 ) 
Frob(L M ) 



k + 1 
4 

k + 2 
6 



h in (z) 

/iim(z) + 10 
k + 2 



k + 1 
3 

k + 



+2 



5 

^22(2) + 



h 2 i (z 

/i 2 n(z) + 4 
fc + 2' 

3 )/l4(Z 



k + 1 
2 

fc + 2 
4 



hi(z) 



Proof of (32). An extension of the "cyclic lemma" of [1] gives a combinatorial proof of 
Identity (32). Consider the action of S n , on pairs (£',{3') in L fc _ ln x P fc _ ln _ 1; that 
permutes labels in £' (as in Section 2). Assume that a' is the underlying path of £'. Since 
£' and /3 are independent, in view of (7), the corresponding Frobenius may be expressed 
as 

(to + k — 2 N 
fc - 1 " 

We are going to show that there is an equivalence relation on the set L fc _ l n x Pj fc _ ljn _ 1 , 
whose equivalence classes are all of size k, each of which containing one and only one 
labelled polyomino. Furthermore, for all pairs of a given class, the set partition associated 
to the labelled path are all the same. Hence, we may conclude that Frob(Lfc jn ) is equal 
to the right hand-side of (35) divided by k, as formulated in (32). 



(35) 
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Figure 8. Only polyomino in the equivalence class of a pair (£, (3). 

To this end, we work with slightly modified versions, denoted by a and j3, of the paths 
a' and (3'. On one hand, a is obtained by simply adding an east step at the end of a'. 
On the other hand, we get /3 by first shifting /3' to the right (sending each vertex (a, b) 
to (a + 1, &)), and then adding a final north step. Hence a now goes from (0, 0) to (k, to), 
while (3 goes from (1, 0) to (k, to). We further denote by £ the labelled path version of £', 
which is directly obtained by carrying over the labels of a' to a. This makes sense since 
the only difference between the two is the perforce unlabeled final step. We underline 
that by construction, £ and /3 end up at the same point, with a final east step for £ and 
a final north step for /3. 

We now consider the two bi-infinite (labelled in one case) paths, denoted by £°° and 
(3°°, respective obtained by repeatedly gluing together copies of the paths £ and (3 at 
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Figure 9. The k transversal intersections of oo-iterates of binomial paths. 



both of their ends. This is illustrated in Figure 9. Two pairs (i[, j3[) and (i^fy), m 
Lfc-i,n x Pfc_ lj7l _i, are considered to be equivalent if they give rise to the same pair of 
bi- infinite paths, i.e. (if , (3f) is equal to (if, f3f) up to a shift of the origin. Equivalent 
pairs correspond to specific portions of the pair of infinite paths. These portions are 
precisely characterized by points (x, y) where the two paths intersect in such a way that 
the i°° reaches (x, y) by an east step, and f3°° by a north step. We then say that i°° and 
have a transversal intersection at (x,y). Equivalent pairs correspond to portions 
of £°° and both lying in the rectangle having south west corner (x — k,y — n) and 
north east corner (x,y), for transversal intersection points p = (x,y). We denote by i p 
and (3 P these portions. 

As shown in [1] the two paths i°° and have exactly k transversal intersections. To 
finish our proof, we need only observe that, for all transversal intersection point p, the 
set partition of i p is equal to that of I. Indeed, the big steps of i p are simply a cyclic 
shift of those of i, together with their labels. □ 

Remark 3. It is interesting to compare Formula (33), evaluated at k = rn, with the 
conjectured expression [4, Formula (17)] for the Frobenius characteristic of the space of 
trivariate generalized diagonal harmonics. Recall that, up to a sign twist, this expression 
takes the form 

53(rn+l)^- a nf (r + 1) 0— ■ < 36 ) 
\hn jex \ J / Zx 

It has been shown in [6] that this formula gives the Frobenius characteristic for the action 
of the symmetric group on labelled intervals in the r-Tamari posets. Our point here is 
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that the difference 

J- (rn)^ ( (r + l) _ n ~ 2 ) ^ - (rn + l)^" 2 ]J ( (r + . 1} j ) ^ (37) 

turns out to be /i-positive. In terms of the corresponding § n -modules, this implies that 
the latter is contained in the former as a "permutation" S n -module. In fact, more appears 
to be true as discussed in Remark 7. 



6. Doubly labelled polyominoes 

We may further add labels to horizontal steps of the lower path, and get an action of 
the product group *§> k x S n on the resulting set ~Lr k n of doubly labelled parallelogram 
polyominoes. We may see that the associated Frobenius characteristic is given by the 
formula 

Frob(L^)(y,z) = -^—h k [(n - l)y] h n [k z] + -L-h k [ny] h n [(k - l)z] 

- (n-l)(k-l) hk[{n ~ 1)y] hn[(k ~ 1)ZL (38) 

Observe the evident symmetry 

Frob(4 2 i)(y,z) = Frob^Xz, y), (39) 

which corresponds to exchanging the coordinates in our polyominoes. For example, we 
have 

Frob(Lg)(y, z) = (6 s 3 (y) + 3 s 21 (y)) s 2 (z) + (3 s 3 (y) + s 2 i(y)) s u (z). 

Here, the s At (y) encode irreducibles of S k , whereas the s\(z) encode those of S n . Naturally, 
our previous formula (32) appears as the multiplicity of the trivial representation of B> k 
in (38). More precisely, recalling that 



(h k [ny},h k (y) ) 



n + k — 1 
k 



we have 

(2) 



Frob(L M )(z) = (Frob(LjgJ(y,z)A(y)) 

h n [k z] 



n — 1 \ k 

fc- 



since the expression in the red parenthesis vanishes. 



1 fn + fe - 2 N 
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It also follows, directly from (38), that the total number of doubly labelled parallelogram 
polyominoes in is 

|Lgj —k n (n- l) fc_1 + (ifc - l)™" 1 n k - (k - l)™" 1 (n - l)"" 1 (n + k — 1). (40) 

(2) 

Considering S^-isotypic components of ^ k n , we sometimes obtain nice formulas for the 
Frobenius characteristic of the resulting S n -module. For instance, denoting by p the 

(2) 

rectangular partition r n , we find that the coefficient of s p (y) in Lrn,n(y, z) is 

M(m-l)z)]. (41) 



li 2 J 



rn.nKy i Z J 

s P [y) rn-1 
See remark 6 for a g-analog of this formula. 

A slight variant holds special meaning in relation with stable configurations for the 
sandpile model associated to the complete bipartite graph K k ^ n [9]. For this variant, we 
fix the label of the first step of the blue path to be 1 (in terms of the sandpile model, 
this step corresponds to the sink). With this convention, the group § n x S^._i acts on 
labels (fixing the label 1 on horizontal steps), and we denote by Lj^ (y, z) the resulting 
Frobenius. 



^;:\y,z) = ^-h n [kz]h k . 1 [ny] (42) 



Proposition 4. We have: 
which implies that 

\L?,: ) \ = k n - 1 n k -\ (43) 

which is the number of K k ^ n -parking functions (or of recurrent configurations, or of span- 
ning trees). 

Proof. We need only add labels 8 to horizontal steps of the blue path in the proof of (32), 
and the argument still holds. □ 



7. Operators and Macdonald polynomials 



To go on with our presentation, we need to recall a few notions related to the theory 
of Macdonald polynomials, and operators for which they are joint eigenfunctions. The 
integral form Macdonald polynomials 9 -£/^(z; q,t), with /i a partition of n, expand 

as 

H^z; q,t) = ^2 K x ,^(q, t) s A (z), 

Ahn 

where the K(q,t)'s are certain rational fractions in q, t, normalized so that K n ^ = 1. 
They are entirely characterized by this normalization together with the two triangularity 
properties 

8 Taking into account the special role of 1. 

9 In previous papers, the notation is often use for these, to distinguish them from other polynomials 
that are not used here. 



LABELLED POLYOMINOES 



17 



(1) H^[z (1 — q)} is in the span of the s\(z), for A >z \i in dominance order, 

(2) H^z (1 — t)} is in the span of the s\(z), for A >z \l' in dominance order. 

It may be worth recalling that we have the special case 



H n (z;q,t)=Y[(l-q i )h n 



i=l 



1-q 



(44) 



Observe that the right hand side is independent of t. This is helpful when we specialize 
the H^'s at t — 1, since in that case we have the following multiplicative property 



H^(z;q, 1) = JJi? Mi (z;g, 1), for /i = /i ly u 2 ■ • • /x r 



(45) 



Thus, we can evaluate this last right hand side using (44). Observe also that this states 
that H^z; q, 1) is proportional to h^[z/(l — q)]. 

Another often used specialization is at t — 1/q. In this case we have (see [15]) 



H,(z;qA/q) = q- n ^ ]J(1 - q 



hook(a) N 



1-q 



(46) 



Operators. For a symmetric function /, we consider the linear operator Af, on sym- 
metric functions, such that 

AfiHj^flBpfaWHp, (47) 

where f\By\ stands for the evaluation of / in the q % t\ with (i, j) varying in the set of cells 
of \x. It is customary to denote by V the operator A £n restricted to symmetric functions 
of degree n. Observe that the operator A hk (resp. A e J appears as the coefficient of u k 
in the linear operator $ such that 

*W = ( II !_L J g M, r esp. *(fl" M ) = ( J] (l + g¥t*))fT M . (48) 

One of the interests of the operator V is that it allows a very compact formulation of the 
Frobenius characteristic of the S n -module of bivariate diagonal hamonic polynomials, in 
the form V(e n (w)) (see [14] for more on this). 

We will also denote Aj (resp. A/) the specialization at t = 1 (resp t = 1/q) of the 
operator Af. It turns out that the operators V, $, and \1/ are multiplicative, since the 
associated eigenvalues are multiplicative. In view of our previous observations (see (44)), 
the set 

{^[ Z /(l - q)] | fi h n} (49) 

constitute a basis of eigenf unctions for the operators Af. Since we are interested in 
calculating the effect of our operators on the elementary symmetric function e n , it is 
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worth recalling that the dual Cauchy formula implies that 



1 — x 



1-q 
1 



1-q 



(50) 
(51) 



where we have denoted by / M (z) the forgotten symmetric functions, that are dual to the 
elementary. We may specialize these expression, setting x = q, to expand the elementary 
symmetric function e n (z) respectively in terms of the s M [z/(l — q)} and /i M [z/(l — q)]. 

Similarly, from (46), it follows that the s M [z/ (1 — g)]'s are eigenfunctions of Ay, since they 
are proportional to the i? M (z; q, l/g)'s. The corresponding eigenvalues are f[B^(q, 1/g)]. 
Observe that in the case of hook shapes, i.e. \i = s^-ryr, these are simply 

f[B {n . r)tl r(q,l/q)] = f[q- r [n} q ]. (52) 

More operators. Recall from [13] that we may introduce symmetric functions E n>r by 
means of the plethystic identity 



1 — x 



E 



{ [q, q)r 



In particular, setting z = q, we find that 

e n (z) = ^E n>r (z). 



(53) 



(54) 



r=l 



These play a key role in a generalized version of the "Shuffle Conjecture", together with 
another family of linear operators on symmetric functions defined as follows. For each 
integer a, an any symmetric function /(z), we set 



C a f(z) 



-qf- a f 



z — 



q-1 

qx 



5> 

k=0 



fc(Z) 33 



(55) 



For example, we have 

C (l) = (-g) 1 ^ s a (z), C ( Sl (z)) = (-g) 1 ^ (s a>1 (z) + q- 1 s a+1 (z)), and 

C a (s 2 (z)) = (-g) 1 ^ ( Sa , 2 (z) + g" 1 s a+1)1 (z)) + g" 2 s a+2 (z)), if a > 2. 

As shown in [], the operators C a 's satisfy the commutativity relations 

g (CfcC'a + C a _iCVf.i) = C a C h + C 6+ iC _i, whenever a > b. (56) 

Recall also that, for a compositions of n into r parts, we may consider the symmetric 
functions recursively defined as 



£ 7 (z) :-- 



C a (E u (z)) if 7 = av, 
if 7 = 0, 



(57) 
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where we denote by the empty composition; and, we have 



E n A Z ) 



(58) 



7gComp(n,r) 



Remark 4. An explicit enumerative description of V(i? 7 (z)) is conjectured, in [13, 
Conj. 4.5], in terms of parking functions whose underlying Dyck has prescribed returns 
to the diagonal. It appears that a similar statement holds for V r (-E 7 (z)), with the 
heights of returns to the diagonal specified by the composition 7. This is what we 
call the r-generalized shuffle conjecture. In the particular case of 7 = n, since 
E n {z) = (—q)~ n+1 h n (z), this conjectures implies that u V r ((— q)~ n+1 h n {z)) is the q- 
graded Frobenius characteristic of the S n -module generated by r-parking functions whose 
underlying r-Dyck never returns to the "diagonal" (of slope except at its end-points. 

Observe that the commutativity relation (56) imply that, for any length r composition 
7 of n, we may express -E 7 (z) as a linear combination of the E^(z)'s, for fi varying in 
the set of partitions of n having r parts. These last symmetric functions are in fact 
proportional to a specialization of Macdonald polynomials, that is 



where, as usual, n(/i) stands fo ^[=i(* — -0 Notice that the right hand side of this 
last identity is expressed in terms of the conjugate partition //. 



8. Formula for the (/-Frobenius of the labelled polyomino module 

Michele D'Adderio has proposed (see [7]) the following explicit formula for the q- weighted 
area counting version of the Frobenius characteristics 



Here A^ fc stands the specialization at t — 1 of the operator A^ fc . This implies that 
Ah k {e n (z)) is e-positive, i.e. with coefficients in N[q\. In generating function format, we 
may reformulate (60) as 



E lx (z) = q- n ^(-l/qr- r H^(z;q,0), 



(59) 



Frob(L fc+ i iri )(z; q) = uA hk (e n (z)). 



(60) 



y^Frob(L fc|W )(z; q) 



u 



k 



en(z))) 



k>l 
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Small values are as follows: 

u 



^Frob(L M )(z;g)M fc = — -h(z), 



k>l 



^Frob(L fci2 )(z;g) 



u 



k>l 



(1 — u)(l — qu) 



$>ob(L M )(z;g) M fc = " 



ho(z) 



1-u 



hn(z] 



k>l 



(1 — u)(l — u q)(l — uq 2 ) 



, . . u (2 + q) , . . u 2 (1 + q) , 



1 - u 



[1-uf 



Even if Formula (60) is still conjectural, the following two formulas should help under- 
stand how uA hk (e n (z)) could be interpreted as a bivariate Frobenius characteristic for 
labelled polyominoes, with a second statistics accounting for the parameter t, on top of 
q which already accounts for area. This also confirms indirectly that (60) should hold. 

Proposition 5. The q-free component of u A/ lfe (e n (z)) is given by the following formula 



id 



A hk (e n (z))\ q=0 = -ih 



l_ t k+i 
1 - t 



- h r 



Moreover, we have 

wA Afc _ 1 (e„(z))| t=1/g 
which is a direct q-analog of (32). 



q 



n+k—nk—l 



1 - q k 
'1-9 



1 - t k 
' 1-t 



n + k-2 
k-1 



(61) 



(62) 



J 9 



For a combinatorial interpretation of the right hand side of (61), see Remark 5. 



Proof. To prove (61), we first observe that A/ lfc (e n (z)) is symmetric in q and t, we may 
thus exchange their role. We may also remove u, so that we need only show that 



A hk (e n (z)) = — [e 



l-q 



k+l 



l-q 



- e r 



i-q 

'l-q 



fel 



(63) 



where Ah k := Ah k \ t _ Q . Evaluating (53) at z — q + , we get 



l-q 



k+l 



l-q 



E 

r=l 



k + r 
r 



(64) 



since, by definition, 



k + r 
r 



_ {q k+1 ;ql 

q (q;q)r 
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Using (64), and well known properties of g-binomial coefficients, we can calculate that 
the right hand side of (63) is 



n 1 

A, fc (e n (z)) = JJ- 



r=l 



r=l 





fc + r 




+ r — f 






r 


9 


r 





^n,r(z) 



A; + r — f 
r — 1 



E n A z ) 



(65) 
(66) 



Now, recalling (54), this last equation readily follows from the surprising fact that the 
_E njr .(z)'s are eigenfunctions of the operator A hk . Indeed, as we will see, we have 

k + r — 1 
r — 1 



A hk E n<r (z) 



E n Jz). 



(67) 



This follows from the even stronger fact that for all composition 7, having r parts, we 
have 

k + r — 1 
r — 1 



A/> fe £ 7 (z) 



E,(z). 



(68) 



In view of the observation following (58), we need only show that (68) holds for E^z), 
with /1 a length r partition. Since we have already seen recalled (see (59)) that these 
E^s are essentially specializations of the H^s, it only remains to check that 

k + r — 1 
r — 1 



A ftfc i2>(z;g,0) 



# M '( z ;<?,o). 



But, by definition of A^,, , this is equivalent to the classical fact that 



k + r — 1 




"1 -g r " 


r — 1 


9 


U-<zJ 



r -lx 



h k (l,q,...,q 



This ends our proof of (61). 

To prove (62), we start by recalling from above (see (46)) that the s M [z/(l — q)]'s are 
eigenfunctions of A/ lfc _ 1 = A /lfc _ 1 \ t _ x . , with eigenvalue hk-i[B^(q, 1/q)]. Once again, we 
remove u, to turn the identity we have to show into 



A fcfc _ 1 (e n (z)) 



9 



n+fc— nfe— 1 



1 - g fe 
; 1-9 



n + k - 2 
fe- 1 



(69) 



As we have already seen, using Cauchy's identity, we have that 



1 -q k 



1-q 



1-q 



Since s v [l — u] is only non zero when v is a hook shape, and S( n - r ) t ir[l — u] = (— u)' r (l — u), 
we may rewrite the above identity clS clS 



1 - q k 

'i-q 



n-l 



£(1 - ?*) (-? fe ) n - r - 1 S( n-0,K 



r=0 



1-q 



(70) 
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Using this identity at k = 1, we may expand e n (z) in terms of eigenf unctions for the 
operator A/ lfc _ 1 , and thus we calculate (using (52)) that 

n-1 



A hfe _ 1 (e n (z)) 



-r-l. 



(n-r),l r 



A/i fc _j £(1 " ?) (-«)" 

r=0 

n-1 

^(1 _ q) (- q y-r-\-^) hk _ x 
n-1 

£(-i) n 



r=0 

n + fc - 2 
fc-1 



-r-l; 



g-1 

1 -q n ~ 
1-q. 



'(n-r),l r 



'(n-r),l r 



q 



n+k—nk—1 



<? r=0 

n + k-2 

k-1 



n-1 



q 



„fc\n-r-l 
"9 J S(n-r),F 



Q r=0 



q-1 



□ 



which shows (69), using (70) to evaluate the sum. 

Remark 5. The right-hand side of (61) may be interpreted as a graded Frobenius of 
ribbon shaped labeled parallelogram polyominoes of size n x k. The needed grading 
corresponds to the dinv statistic for unlabeled parallelogram polyominoes, as defined in 
[2]. Recall that, in the A- word encoding w = u>iu>2 • • -w/- of a parallelogram polyomino 
(see Figure 5), the dinv statistic is the number of pairs (wi,Wj), with i < j and j is the 
successor of % in A. For any ribbon-shaped labeled polyomino, we may simply set its 
dinv-statistic to be equal to that of the underlying parallelogram polyomino. With this 
definition, we observe that the dinv-value coincides with the area below the lower path 
of the polyomino. This readily implies that 



Frob(L fcin )(z; q) 



t n h 



1-t 



h, 



1-t 



(71) 



Remark 6. Experiments suggest that we have the following area enumerating analog of 
Formula (41) 



s p(y) 



uV r {{-q)- n+1 h n 



(72) 



where p is the partition consisting of n rows of size r. In fact, this identity follows 
from the r-generalized shuffle conjecture (see Remark 4), considering paths that never 
come back to the diagonal of slope 1/r. It would be interesting to describe all isotypic 
components of Lrn,n in this manner. 



Toward a general dinv statistic. A dinv-statistic 10 on general labelled polyominoes 
would make possible an analog of the "Shuffle Conjecture" of [12], in the form: 

A hk (e n (z)) = J2q aTCai7T) t diQv{n) z^ (73) 

7T 

with 7r varying in the set ~hk,n- More explicitly, this means that we would consider more 
general labellings of polyominos; with labels in N and allowing for repetition of labels, 



Which we don't know how to define yet. 
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while keeping the strictly increasing condition alogs runs of vertical steps. In the above 
expression, z n denotes the product of the variables Zj, for i running through all labels 
(with multiplicities) of n. In a sense, this should be a natural extension of the original 
Shuffle Conjecture, since it appears (observed experimentally) that 



A ftrn _ x (e n (z)) - )V r (e n (z)) lies in N[g,t]{s„|/J h n}. 



(74) 



As we have already seen (see Remark 2), parking functions can be considered as special 
cases of labelled parallelogram polymominoes, making (74) a natural phenomenon. 

Remark 7. Explicit computer calculations suggest that even more seems to hold. In- 
deed, up to the same multiplicative factor of q\ 2 ), the difference between A hn _ 1 (e n (z)) 
and the "Shuffle-like Conjecture Formula" for the Frobenius of the space of trivariate 
diagonal harmonics (see [1, Formula (15)]), appears to always be Schur-positive (with 
coefficients in N[q, t]). Not only is this stronger than (74), but it suggest that a third 
statistic might play a role here. 



Remark 8. Along these lines, it is interesting to observe that, when k < n, 



A P Je n (z)) 



1 



k 



k+l 



1 



Hence, the Frobenius (10) may be expressed as 

Frob(L M )(z;g)=g fc "-(^ 1 )[A;+l] 
Observe the similarity between this expression and (60) 



q-1 



uA ek (e n (z)). 



(75) 



(76) 



The bounce statistic. Angela Hicks has recently given a nice bijection from which it 
follows that we have the following explicit formula due to Garsia and D'Adderio 



P*,»(g,*) := 



= (V{e k+n „ 2 ),h h -iK-i) 



(77) 



where b(n) denotes the bounce statistic of a polyomino tt, introduced in [9]. In other 
words, up to multiplicative factor, the right hand side is the coefficient of the monomial 
symmetric function of m/-_i in _i in the expansion of V(efc +n _2). It immediately follows 
that Pfc, n (<7, t) is symmetric, both in q,t and n, k, as already observed in [9]. Moreover, 
specializing t at 1/q, we get the following direct g-analogs of (17) 



q (k-l) (n-X) p ( ^ l/q) 



[n + k] q 



n + k 




~n + k - 2" 




n 




n — 1 






g 




9 


;-l,fc-l)( 


1. 


q,q 2 ,...,q n 


)• 



(78) 
(79) 



To prove this, we proceed as follows. Recall from [15] that 

1 



t<-i/q [n + 1 



-e n [z [n+ 1 
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so that (77) may be reformulated as 

7 — — r — jY e n+k-2 [z [n + k - 1 ) g ) , h k -i h 



/ n+ fc_2\ / 1 

gV 2 JP M (g,l/g) 



n— 1 



Now, Cauchy formula gives 

e n+fc _ 2 [z[n + A;- 1],] = s A ([n + fc - 1 ],) s A (z) 

Ahn+fc-2 

leading us to 

1 



qy * )¥ kin {q,l/q) 



n + k — 1 



-ejfe_i([n + fc - l] q ) e n -i([n + k - l] q ). 



Finally, since 



we conclude that 

(n + k-2\ 

gl ^ )¥ k>n (q,l/q) 



e A\ n ] q ) = (f 



n + k — 1 



-g 



fe-i 



n + k — 1 




n + k — 1 


k-l 


9 


n — 1 



JO) 



51) 



(82) 



A direct calculation shows that this is equivalent to (78). 

9. Thanks and future considerations 

We would like to thank Vic Reiner for reminding us that we should consider doubly 
labelled polyominoes in the specific manner we describe here. We would also like to thank 
Nolan Wallach for putting us on the right track for the interpretation of parrallelogram 
polyominoes as indexing sets for .SZ^-invariants. 

In our view, the main question that begs to be addressed in the future, seems to find a 
natural context in which A hk (e n ) would appear as the bigraded Frobenius characteristic 
of an § n -module similar to the S n -module of diagonal harmonic polynomials. This would 
help truly explain all our observations above. 
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